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For the antiferromagnetic J\-Ji quantum spin chain with an even number of sites, the point 
4 = Ji/2 is a disorder point. It marks the onset of incommensurate real space correlations for 
J2 > J2 • At a distinct larger value of J2 = 0.52036(6) Ji , the Lifshitz point, the peak in the static 
structure factor begins to move away from k = ir. Here, we focus on chains with an odd number of 
sites. In this case the disorder point is also at Ji = Ji/2 but the behavior close to the Lifshitz point, 
J2 ~ 0.538Ji, is quite different: starting at J2 , the ground state goes through a sequence of level 
crossings as its momentum changes away from k — tt/2. An even length chain, on the other hand, is 
gapped for any J2 > 0.24Ji and has the ground state momentum k = 0. This gradual change in 
the ground state wave function for chains with an odd number of sites is reflected in a dramatic 
manner directly in the ground state on-site magnetization as well as in the bi-partite von Neumann 
entanglement entropy. Our results are based on DMRG calculations and variational calculations 
performed in a restricted Hilbert space defined in the valence bond picture. In the vicinity of the 
point J2 = Ji/2, we expect the variational results to be very precise. 



I. INTRODUCTION 

Disorder points were first discussed by Stephenson in 
models described by classical statistical mechanicsP^ On 
one site of a disorder point, the correlation function shows 
monotonic decay, on the other oscillatory decay. Depend- 
ing on the how the wavelength of the oscillation depends 
on the temperature, one distinguishes between two kinds 
of disorder points. If the wavelength of the oscillation 
depends on the temperature, one speaks of a disorder 
point of the first kind, if it does not, one speaks of a 
disorder point of the second kind. 2 ' In the first studies, 
disorder points were found where the paramagnetic phase 
of frustrated two-dimensional Ising models starts to show 
incommensurate instead of commensurate behavior. In 
models with competing commensurate and incommensu- 
rate order one might expect such a point to occur where 
the short-range correlations with the largest correlation 
length change from being commensurate to being incom- 
mensurate. Such a point should then be associated with 
a cusp in the correlation length, a fact that was quickly 
established. 5 Schollwock, Jolicoeur and Garel first inves- 
tigated disorder points in a quantum spin chain for the 
bilinear-biquadratic S = 1 quantum spin chainj^l which 
has H = cos 9 Si • Si+i +sin 6 Si-Sj+i. They pointed 
out that the disorder point in this gapped quantum model 
coincides with the AKLT-point where tan 9vbs = 1/3 
and the known ground state is a valence bond solid (VBS) 
with correlation length £ = l/ln(3). They also identified 
another distinct point, the Lifshitz point, at tan#^ = 1/2, 
where the peak in the structure factor is displaced away 
from 7r due to incommensurability effects. A third distinct 
point in this S = 1 model, taitOdisp — 0.4, has also been 
located^ where the minimum in the magnon dispersion 
shifts away from tt and the curvature (velocity) vanishes. 
These 3 points can be distinct since no phase transition 
occurs and the correlation length remains finite. Subse- 
quently it was confirmed 8 that for the S — 1/2 J\-Ji spin 



chain with Hamiltonian: 
H = Ji s * • s *+i + J 2 S * • S «+2> J i> J 2 > . C 1 ) 

i i 

the situation is similar. For the calculations presented 
in this paper, we set J\ = 1 and vary the remain- 
ing parameter Ji. The disorder point, with minimal 
correlation length (£ ~ 0), occurs at the Majumdar- 
Ghosh 9 (MG) point, J^ — Ji/2, and the Lifshitz point at 
J 2 L = 0.52036(6) JiP 
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FIG. 1. (Color online.) The odd length J1-J2 chain. Two 
different dimerization patterns are separated by a soliton. 

At the disorder points, the ground states of these two 
quantum spin models share important features: the sys- 
tem is gapped and an exact wave function is known. For 
both, the momentum of the lowest excitations changes 
at a distinct point. Yet, there are also important differ- 
ences between the two systems. While the S = 1 VBS 
state remains an exact state for a chain with an odd num- 
ber of spins, this is not the case for the 5=1/2 J\-J<i 
chain at the MG point where no analytical expression 
for the odd length ground state wave function is known. 
Moreover, the odd length J\-Ji chain is gapless in the 
thermodynamic limit within the Sxot = 1/2 subspace, 
and a large spin gap exists. The onset of incommensu- 
rability effects for odd length chains must then be quite 
distinct from the onset in even length chains. Here, we 
show that this is indeed the cases. While the disorder 
point remains unchanged, the nature of the Lifshitz point, 
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At J.f, a sequence of 



2 — 0.538Ji, is rather different 
level crossings starts, changing the ground state momen- 
tum away from k = tt/2. Although the correlation length 
remains small, the change in ground state momentum 
induces pronounced oscillations directly in the on-site 
magnetization as well as the entanglement entropy. The 
modulations in the on-site magnetization are potentially 
observable in experiments. This scenario is reminiscent of 
a real Lifshitz transitioii^L 3 -' in which the ground state 
becomes modulated. The scaling of the entanglement en- 
tropy at L ifshitz transitions recently has been the subject 
of interest)!^ 

The S — 1/2 J1-J2 antiferromagnetic (AF) spin chain 
is one of the simplest frustrated Heisenberg spin models, 
but it has a rich phase diagram. The system undergoes a 
transition^ from a gapless Luttinger liquid to a dimer- 
ized phase at a critical value of Jf = 0. 241167 Ji^^l 
For even length chains the ground state wave function 
at the MG point J2 = J\J2 is known to be formed by 
nearest neighbor dimers P 20 * 21 ! It is two-fold degenerate, 
corresponding to the two possible nearest neighbor dimer- 
ization patterns, indicated in Fig. IT] A s is evident from 
Fig. [l] an unpaired spin, a solitonp^2l can a ct as a 'do- 
main wall' and separate regions of different dimerization 
patterns. In the Luttinger liquid phase unpaired spins are 
more commonly called spinons since they do not act as 
domain walls. Spin excitations in the even length chain 
correspond to introducing two solitons and it is knowrPS 
that in the vicinity of the MG point the solitons do not 
bind and a large spin gap of A = 2A so / (at the MG 
point A soi /Ji = 0. 1170(2) ]P exists. The spin gap for 
even length chains is known to remain sizable^ES] beyond 
J%. The presence of a large soliton mass, A so i, renders 
variational calculations based on a redu ced H ilbert space 
consisting of soliton states very precise J22122I a fact that 
we shall exploit here. 

In contrast, for odd length chains it is not possible for 
the chain to be fully dimerized and the ground state wave 
function is not known for any value of J2- An 5=1/2 
soliton that effectively behaves as a free particle^ is 
always present in the ground state and gives rise to gapless 
excitations. Depending on the quantity in question, odd 
and even length chains can show very different behavior. 
Under open boundary conditions (OBC), this has for 
example already been seen in th e on-site magnetization,^ 
the entanglement entropy^^U and the negativity!^] As 
mentioned, here we focus on odd length chains. 

While it is possible to perform highly precise density ma- 
trix renormalization group (DMRG) calculations we ll be- 
yond the onset of incommensurability for even chains J25J2SI 
the sequence of level crossings that we encounter for odd 
length chains for J2 > significantly restrains the use- 
fulness of the DMRG technique in a large region of pa- 
rameter space for J 2 > J% ■ Fortunately, using the picture 
of Shastry and Sutherland,^ it is possible to quite effi- 
ciently perform very precise variational calculations for 
both open and periodic boundary conditions (PBC). Here, 
we mainly present results of such variational calculations 



and supplement them with DMRG-results. 

A number of spin-Peierls compounds, which to some 
extent realize the J\-Ji spin chain, have been identified. 
One of the most well known is CuGeOaP^In these ma- 
terials impurities often cut the chains at random points. 
Therefore both odd and even length chains are present. 
A particular point of focus has been the study of S = 1/2 
solitonsSSHSSl j n these systems. Thus, our results might 
be directly verifiable if materials with sufficiently large 
J 2 > J 2 can be found. 

The outline of the paper is as follows. In section [TT| the 



variational approach is described. Section III begins with 
a presentation of our DMRG results for the correlation 
functions, correlation lengths and the structure factor. In 
section [Til A| we discuss our variational results for the J\- 
J2 with periodic boundary conditions and show the change 
in ground state momentum developing at the Lifshitz 
point. Section |III B| contains variational and DMRG 
results for the on-site magnetization and level crossings 
occurring with open boundary conditions. Variational 
and DMRG results for the entanglement entropy for a 
range of J2 for odd length chains (OBC) are presented 
in section |IV| and contrasted with results for even length 
chains (OBC). Finally, estimates for the location of the 
Lifshitz point are presented in section |V} 

In the following, we shall take Ji = 1. This leaves us 
with only one parameter, J2, that governs the properties 
of the system. 



II. THE VARIATIONAL METHOD 

Most of the results presented in t his paper w ere gen- 
erated using variational calculations) 22 1 23128 1 30 1 36 1 i.e. the 
results were obtained by minimizing the expectation value 
of the Hamiltonian within a reduced Hilbert-space: 
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and the minimization is done with respect to the Cj . To get 
a good estimate of the true ground state of the system, it 
is necessary that the ground state has a sizable projection 
onto the subspace one diagonalizes in. The quality of 
the result of a variational calculation thus depends very 
strongly on the choice of subspace. Often one has to 
rely on physical insight and intuition to choose well. For 
the J\- J2-model, which we consider, the selection of an 
appropriate subspace is straight-forward as long as one 
stays in the dimerized phase. In contrast, in the Luttinger 
liquid phase, selecting an appropriate subspace seems 
intractable. 

The first variational calculations on the J\- J2-model 
we re don e in a space that we in the following shall call 
Jl 1 22 | 23 | jf. j g gpanned by the states in which there are 
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domains that have one of the two ground state configu- 
rations of the MG-chain and which are separated by one 
soliton. Examples can be seen in Figs. [T] and [2] The 
arrows in Fig. [2] serve to fix the phase of the dimers that 
make up the ground state. Our convention is such that 
if the arrow goes from site i to site j the spins are in 
the state: 4g(| i)j - I t)i)- F° r a chain with 
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FIG. 2. Two variational states used in the calculations within 
Ro for a chain with an odd number of sites. The arrows 
between sites are used to fix the phase of the dimers (see main 
text) . 

an odd number of sites, a set of single soliton states can 
be generated by leaving the chain maximally dimerized 
and taking the remaining site to be in the S' z =l/2-state. 
For a chain with open boundary-conditions, the soliton 
can only reside on every second site. The dimension of 
this variational subspace is then D=(N + l)/2. We use 
N to denote the length of the chain. For calculations on 
odd length chains with periodic boundary conditions it 
is necessary to allow a nearest neighbor dimer across the 
boundary and to let the soliton cross the boundary by go- 
ing from site N to site 2. In this case Ro has dimension N 
and incorporates states with the soliton at every site with 
the remaining spins paired in nearest neighbor dimers. 
(For odd N and PBC it becomes difficult to distinguish 
the 2 dimerization patterns since they twist into each 
other at the boundary. Still, the soliton clearly denotes a 
'domain wall' between the two patterns). 

To improve upon Ro, it is natural to act with the Hamil- 
tonian onto the space as doing this repeatedly generates 
a space that contains the ground state if the starting 
space had any overlap with and all symmetries of the 
ground state. It was shown that acting onto Rq with 
the Hamiltonian only once is at the MG point already 
enough to make the calculation almost exact For the 
Ji-J2-model with Ji ^ 0.5Ji the linearly independent 
states generated by acting with the Hamiltonian onto 
Ro fall into three classes, each of which corresponds to a 
variational subspace: 

• The variational space R : Spanned by the states 
that are in Rq. 

• The variational space R\: Spanned by states in 
which sites to the left and right of the soliton are con- 
nected by a valence bond. Pictorial representations 
of example-states are shown in Fig. [3] 

• The variational space R2: Spanned by states in 
which two neighboring sites are in a valence bond 
with their next-nearest neighbor. These states are 



generated by the action of the nearest-neighbor- 
terms and the next-nearest neighbor terms in the 
Hamiltonian on adjacent dimers in the states in 
Rq. Pictorial representations of example-states are 
shown in Fig. |4j In the case of the MG-chain, J\ 
and J2 are balanced in such a way that these states 
are not generated because they occur with a weight 
of 2J 2 - Jr. 

The number of states in Ro and R\ scales linearly with 
the size of the chain, whereas the number of states in 
R2 scales quadratically. Due to computational cost we 
have thus not found it practical to use the union of the 
three as the variational subspace for chains longer than 
101 sites. We performed calculations using the union of 
i?o and i?i (in the following called Z s ) for chains up to 
1001 and the union of all three (in the following called 
Zb) for a chain of 101 sites. In this way we could go 
to long chains and also check the validity through the 
comparison at N = 101. We found that while there were 
small quantitative differences between calculations done 
in Z s and Z^ the overall qualitative features of the results 
where the same. Therefore, we chose to use Z s , the union 
of Ro and R x , or just R for the variational calculations 
shown in this paper. 

All the states in these spaces have S^ ot := J^i &l = 1/2- 
We could equally well have worked in the S^ ot = — 1/2 
space. States of higher total spin are of little importance 
to the low-energy physics since they contain more solitons 
and are thus gapped by at least twice the soliton mass 
A so ;. Since A so / is sizable^ in the regime of our study, 
such states can be disregarded for both odd and even N. 




(b) 

FIG. 3. Two variational states used in the calculations within 
i?i for a chain with an odd number of sites. The arrows 
between sites are used to fix the phase of the dimers (see main 
text). 

A variational description of a chain with an even number 
of sites can be done along the same lines. Again, the states 
are chosen in order to leave all but two spins in the favored 
dimerized state. In this way one can gain insight into 
the low-energy singlet- as well as triplet-excitations by 
choosing the two spins to be in singlet- or the triplet-states 
respectively. 

If one considers a subspace with an orthogonal basis, 
one can just diagonalize the Hamiltonian. While an easy 
way to orthogonalize Rq is known,^ this is generally not 
true for other subspaces. Importantly, for Ri no such 
method is known. We thus have to solve the generalized 
eigenvalue-problem given by 

Si = A 580 , (4) 
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FIG. 4. Three states that are generated by acting with the 
Hamiltonian onto Ro- The arrows between sites are used to 
fix the phase of the dimers (see main text) . 



where Sjij := (ifi\Htpj) and 2$ij := ((pi\tpj). Such 
generalized-eigenvalue-problems can be solved numeri- 
cally by standard routines. We calculate f) and 03 by 
evaluating their defining expressions. This is possible 
because for valence-bond-states the action of H on them 
as well as the overlap between them can straightforwardly 
be calculated in an automated manner. How to do all 
other calculations necessary to get the results presented 
in this paper has already been described in an earlier 
publication.^! We took the coefficients Cj in Eq. ^ to be 
real. Also for PBC, the resulting wave function is not an 
eigenstate of the translational operator which would have 
required the use of complex Cj's. Effectively we obtain 
states that are linear combinations of translationally in- 
variant states with k and —k, degenerate in energy. While 
this has no bearing on the obtained energies, it affects 
real-space quantities like the on-site magnetization and 
entanglement which cannot be translational invariant. 



O-10 



o 



-20 



(a) 



I I 
o 0.5 

o 0.51 



0°o 



5 
4 

e 3 
2 
1 



.(b)' 



in 



Odd 
Even 



So 





10 20 30 
i 


0.495 0.510 

J 2 




(C) Odd ^ 


l . — i — | l l 
| Even| 


1.2 


// — 


// ~~ ~ 


0.6 


// 


AT 




0.51 


0.51 


0.0 


0.57" 

i i i i 


0.57" 

i i i i 



0.2 0.4 0.6 0.8 

q/ix 



0.2 0.4 0.6 0.8 

q/ix 



FIG. 5. (Color online.) (a) The spin-spin correlation function 
C(r) — {SiSi+ r ) for a chain with 201 sites and J2 = 0.5,0.51. 
The correlation-length (b) and the static structure factor (c) 
for chains with an odd (201) or an even (200) number of 
sites. The data was obtained with DMRG. Open boundary 
conditions were employed and m = 256 states kept. For both, 
odd and even length chains, the correlation-length displays 
a minimum at the MG point. The static structure factors of 
odd and even are shown for J2 = 0.51 and J2 = 0.57. The 
maximum remains at q = n until the Lifshitz point is crossed 
(not shown in figure). 



III. THE INCOMMENSURATE BEHAVIOR 

Previous num erical stu dies of disorder points in S = 

!/ gnmrnmim and s = 1 mmm quan t um SP m chains 

have concentrated on the behavior of even length chains. 
For the S — 1/2 J\-Ji chain it has been shown that the 
disorder point of this 1-dimensional quantum system can 
be understood as a 1 + 1 dimensional classical disorder 
point. In particular, it was 

showrPHUl 

that in the "com- 
mensurate" region of the phase diagram the correlation 
function behaves asymptotically, with r 3> £, as 

(SiSi+r) ~ (-1)^-^, (5) 

and in the "incommensurate" region of the phase diagram 
as 

— r/4 

(SA+r) ~ ( l) r — cos [(«? - n)r + 0] . (6) 
y/r 

Here, q is the wave vector of the incommensurate correla- 
tions and 4> a phase shift. However, right at the disorder 



point separating commensurate and incommensurate cor- 
relations the correlation function is asymptotically purely 
exponential: 

(SA+r) ~ (-lfe-^. (7) 

For these quantum spin models it appears that this purely 
exponential behavior is in part connected to the fact that 
the ground state is an exact nearest neighbor dimer state. 
Interestingly, as we shall see, the correlation functions 
at the MG point for odd length chains display the same 
behavior in the absence of a unique near est ne ighbor 
dimer ground state. Furthermore, it is knowrP2HSl that as 
the disorder point is approached from the commensurate 
side, the derivative of the correlation length with respect 
to the driving coupling becomes infinite, while it is finite 
on the incommensurate side. It is also known that the 
disorder point has special degeneracies that are exact for 
any system size N. For instance, for the J1-J2 chain with 
periodic boundary conditions and an even length, the two 
dimerization patterns are degenerate at the disorder point 
while their symmetric and antisymmetric combinations 
are split with an exponentially small gap away from this 
point. 
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We now present our results for the incommensurate 
effects in odd length S — 1/2 Ji-</ 2 chains. Our first point 
of focus is the location of the disorder point. As stressed 
above, when TV is odd, the nearest neighbor dimer wave 
function is not an exact solution 22 and there are also 
no special degeneracies. There is therefore no reason to 
expect that the behavior of the correlation length at the 
MG point is in any way unique. However, as we shall see, 
this is indeed the case. DMRG results for C(r) = (SiSi +r ) 
for an open chain with N = 201 are shown in Fig. |5ja) 
for J 2 = 0.5 and 0.51. The correlation function follows a 
purely exponential decay at the MG point, J 2 = 0.5, with 
a finite correlation length: 

£mg~2.8 (odd iV). (8) 

Distant spins in odd chains are correlated even at the 
MG point, because the soliton is present in the chain. 
The correlations can be thought of as correlations in 
the soliton wave function. Secondly, as can be seen in 
Fig. |5ja), incommensurate correlations are clearly present 
for J 2 — 0.51. They were present in every calculation 
we performed with J 2 > 0.5. We conclude that the 
disorder point remains at J 2 = Ji/2 albeit with a finite 
correlation length compared to the case of even N where 
the correlation length is nominally zero. 

The precise behavior of the correlation length around 
the disorder point J2 = 1/2 appears to have been studied 
neither for even length nor for odd length chains. Results 
for larger J 2 > 0.6 are available for even iVPBy fitting 
DMRG results for chains of 200 and 201 sites to the forms 
Eqs. ([5]), and Q we have determined £ as a function of 
J 2 for both even and odd N (see Fig. |5jb)). The results 
for the even and the odd length chain are remarkably 
similar. At the disorder point there is a discontinuity in 
the slope of £ and on the commensurate side the slope of 
£ approaches —00. We found that close to the disorder 
point in the commensurate region combined forms like 
\C(r)\ ~ Cexp(-r/^ c )/V^ + Dexp(-r/^ D ) with £ c > 
£d fit the data better than the single forms Eq. §5j, 
and (J6j> , because the dominant short-ranged correlations 
change at the disorder point. The results presented in 
Fig. [5] do not use such combined forms. We also note 
that for odd N and for a range of J 2 > 0.538 it becomes 
very difficult to obtain reliable DMRG results due to the 
appearance of many almost degenerate states. 

The structure facto r for even chains has been studied in 
some detail previouslyPSl anc ; the Lifshitz point has been 
located, J 2 L = 0.52036(6). 8: Our DMRG results are shown 
in Fig. [5^c) . In agreement with previous studies for even 
N, we observe that the maximum in the structure factor 
remains at q = ir for J 2 = 0.51 but has clearly moved away 
from 7r at J 2 = 0.57. This is clearly also the case for odd 
N. Due to the above mentioned difficulties in obtaining 
reliable DMRG results for odd N and J 2 > 0.538 we have 
not been able to determine the precise point where the 
peak in the structure factor is displaced from q — ir. Using 
the variational techniques outlined above it is possible to 
understand in detail what happens close to J 2 ~ 0.538. 
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FIG. 6. (Color online.) Comparison of the estimate by Shastry 
and Sutherland^!, the variational method and exact diago- 
nalization (ED) dataP^ Data was taken for a chain with 23 
sites and J2=0.5. Only the lowest energy of a spin 1/2 exci- 
tation for every momentum is shown. The deviations to ED 
occur at higher energies where the dispersive mode enters the 
continuum and the variational calculation is only of limited 
value. 

A. Variational Results in Periodic Boundary 
Conditions 

We now turn to a discussion of our variational results 
obtained using the method outlined in section |TT] We 
begin by focusing on the case of odd length chains and 
periodic boundary conditions. The case of open boundary 
conditions will be the subject of the next subsection. The 
results shown in this subsection were obtained using the 
space i?o (see section |TT]) , consisting of all single soliton 
states with S*J ot = 1/2. 

At the MG point the spectrum of the Ji-J 2 -model 
has been studied extensively. The feature that is most 
important to us is the low-lying dispersive line that is 
well separated from the continuunPH and roughl y follo ws 
a cosine as found in previous variational studies J22J45I 

E{k) = \ (5 + 4 cos(2fc) - 37V) . (9) 
8 

Our variational method reproduces this estimate and 
agrees well with the low-energy data of an exact diago- 
nalization of a chain of 23 sites (see Fig. [6]). It may be 
surprising that the minimum of the dispersion relation 
is not at k = it but at k = tt/2. This is a natural con- 
sequence of the effective doubling of the unit cell that 
occurs because the action of the Hamiltonian displaces 
the soliton by two sites. 

One of the strengths of the variational method is that 
within the limits of the approximation it is possible to 
easily access not only the ground state but also the en- 
tire energy spectrum within the variational subspace of 
S^ ot = 1/2 states. Computing the spectrum through 
the transition region reveals very surprising behavior (see 
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FIG. 7. (Color online.) The excitation spectrum with periodic 
boundary conditions in the variational subspace with Sj ot = 
1/2. The energies of the first few excited states are shown. 
The ground state energy was set to zero. The energies of the 
excited states approaches the ground state energy until, at 
J2 = J 2 , level-crossings start to occur. All energy- levels are 
doubly degenerate. The data was taken for a chain with 401 
sites. 



Fig. All the states are two- fold degenerate corre- 
sponding to the energetically degenerate k and —k. As 
one approaches the transition, the excited states linearly 
move closer and closer to the ground state. At the Lif- 
shitz point J 2 ~ 0.53 the energy of the first excited state 
crosses the ground state energy. This level crossing marks 
the first shift in the ground state momentum and is fol- 
lowed by a series of other level crossings at larger J 2 that 
further shift the ground state momentum. Clearly, the 
presence of the many adjacent level crossings hinders the 
effectiveness of DMRG calculations. 

This is in stark contrast to the spectrum of even length 
chains: The ground state of even length chains is exactly 
two fold degenerate at the MG point for any N whereas for 
larger J 2 the symmetric and antisymmetric combinations 
are split with an exponentially small gap in N. The 
excited states are separated from these two states by 
a large gap of approximately 2A so [. This gap persists 
throughout the tr ansition region and no level crossings 
are observed! 17 * 26 ! 

It is very instructive to look at how the dispersion rela- 
tion in Fig. Revolves with J 2 . As can be seen in Fig.|8j the 
dispersion relation changes its shape when J 2 is increased. 
The minimum at k = tt/2 first becomes flat very close to 
J 2 = 0.53 and then becomes a local maximum. In the 
process two minima are created, which move away from 
k = tt/2 with increasing J 2 . The ground state momen- 
tum is then clearly changing away from k = tt/2 beyond 
J 2 = 0.53 and we may identify the poin t where this hap- 
pens with a real Lifshitz transitio n 1 11 ! 46 * ^ as opposed to 
the corresponding point in the S = 1 bilinear biquadratic 



FIG. 8. (Color online.) Dispersion relation for varying J2 . The 
data was taken for a chain with 401 sites. To avoid cluttering, 
the dispersion relations for the smaller two J2 were shifted. 
As J 2 is increased, the minimum first flattens and then turns 
into a local maximum. 



chain where the ground state momentum remains un- 
changed and the shift is in the excited magnon dispersion. 
Due to the shift in the ground state momentum we con- 
clude that the maximum of the structure factor will shift 
away from k = tt at the same point. This is consistent 
with the data in Fig. [5] We therefore in the following 
refer to this point as the Lifshitz point, J 2 '. The precise 
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FIG. 9. (Color online.) The second-order coefficient of the 
dispersion relation c changes sign ((a), (b)) while the forth- 
order coefficient d stays positive (b) . The data are results off a 
fit to data of the kind shown in Fig. [8] The dotted horizontal 
line in (a) separates positive from negative values. 

behavior of the dispersion relation close to J 2 is analyzed 
in Fig. [| 

For a range of J 2 , we fitted the dispersion relation to 
the forrrP 

E(k) = E(k ) + - fc ) 2 + ^(k~ fco) 4 (10) 
and confirmed that the second-order coefficient c changes 
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its sign at a J2 close to 0.53 while the fourth order coeffi- 
cient d stays positive (see Fig. [9| . This behavior is typical 
of a Lifshitz transition and if the coefficient c = v 2 /A so [ 
is associated with a velocity, v, the Lifshitz point signals 
the vanishing of this velocity.^ 

The variational calculations with periodic boundary 
conditions presented in this section were limited to the 
subspace Rq described in section |TT] This basis only in- 
cludes nearest neighbor valence bonds and it is quite 
noteworthy that the physics of the Lifshitz point along 
with the associated level crossings are captured within 
this simple basis set. However, as we discuss in section [V] 
we do not expect the precise location of the Lifshitz point 
to be accurately determined within i? - 



B. Open boundary conditions 

In materials that realize the J1-J2 spin chain impuri- 
ties are always present. They often act as non-magnetic 
impurities effectively breaking the linear chains into fi- 
nite segments. The use of open boundary conditions is 
therefore closer to the experimental situation than the 
use of periodic boundary conditions. Furthermore, it is 
natural to expect half of the chain segments to have an 
odd number of sites. In this subsection we therefore focus 
on odd length chains with open boundary conditions. In 
particular, we describe the change that switching from 
periodic to open boundary conditions causes. The vari- 
ational results shown in this subsection were obtained 
using the space Z s (see section [TTj) . 

One quantity that very directly shows the qualitative 
difference between PBC and OBC is the ground state 
on-site or local magnetization, (Sf), whi ch is of impor- 
tance to for instance NMR measurements!^!!] Figure 10 
shows the on-site magnetization at 8 different values of 
the frustrating interaction J2 between J2 = and J2 = 00 



in a chain of 601 sites. The figures 10 c)-(f) show varia- 
tional calculations through the Lifshitz point J% where 
DMRG calculations are less effective, the remaining re- 
sults (a),(b),(g) and (h) are obtained with DMRG. 

The Luttinger liquid phase (J2 < J%): The transition to 
the dimerized phase occurs cit i/^j 

see Fig. |10[b). At this 



point, as well as throughout the Luttinger liquid phase 
(J2 < J£)i the on-site magnetization agrees very well with 
the prediction for the on-site magnetization in the ground 
state with Sj ot = 1/2 from conformal field theory!^ 



(S?) = C(- 




1 

27V' 



(11) 



(St) 



0.03 
0.00 

-0.03 
0.03 
0.00 

-0.03 
0.03 
0.00 

-0.03 
0.03 
0.00 

-0.03 
0.03 
0.00 

-0.03 
0.03 
0.00 

-0.03 
0.03 
0.00 

-0.03 
0.03 
0.00 





-0.03 - 



I h 



H h 



•(c) 



H 1 1 1 h 



■(d) 



H 1 1 1 h 




•(f) 



+ 



+ 



+ 



+ 



+ 



j 2 = 0.535 

.ooo<X>0060060<>o<*>~ 



(g) 



+ 



+ 



+ 



+ 



= 1 .000 



.(h) 



Jo = OO 



J l_ 



_l_ 



100 200 300 400 500 600 
i 



FIG. 10. (Color online.) The on-site magnetization at a range 
of values of J2 for TV = 601 sites. Additional structure appears 
beyond jj". For J2 = 0, Jf, 1 and 00 DMRG-data (obtained 
with m = 256 states kept) are shown (red). For the remaining 
J2 variational calculations are shown (blue). 



where C is a constant. In this phase (Sf) increases with 
the characteristic behavior (Sf) ~ vi for small i close to 
the boundary. 

Dimerized phase with J2 < Jjf' Once the dimerized 
phase is entered (Sf) is drastically altered. The on-site 
magnetization roughly follows the behavior of a massive 



particle in a boxP^I with (Sf) ~ i 2 close to the bound- 
ary. This behavior is clearly visible at the MG point 
(Fig. 10 'c)). As J2 is increased beyond the MG point 
towards the Lifshitz point J-f the central peak sharpens 
(Fig-pTd)). 



'Incommensurate ' phase J2 > J 2 ■ At the Lifshitz point 
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there is another dramatic change in (Sf): additional 
maxima develop and the magnetization is modulated by 
an oscillating function (Fig. [lO^e)). Upon increasing J 2 
further, more such maxima form and the wave-length of 
the modulation decreases (see figures 10 'f) and (g)). If 
J2 is fine tuned for a given N it is possible to find a point 
where 2 maxima occur in (S*), then 3 maxima and so 
forth. 

It is natural to expect this behavior based on the results 
for PBC presented in section |III A| The local magnetiza- 
tion is effectively modulated with the momentum of the 
ground state. The running wave found under periodic 
boundary conditions is converted to a standing wave un- 
der open boundary conditions. Then, as the momentum 
of the ground states changes with growing J2, the wave- 
length of the modulation shrinks. Finally, in Fig. [TO^Ii) 
we show results for J2 — > 00. In this limit the odd length 
chain with N sites is split into 2 chains with (N — l)/2 
and (N + l)/2 sites one of which will have an even number 
of sites and hence (S?) = 0. The on-site magnetization 
of the other chain can be found by calculating (Sf) for 
a chain with J2 = of the same length. The results 
shown in Fig. [l0|h) were obtained in this way, i.e. from 
data for a chain with N = 301 and J2 = that was then 
interspersed with zeros from the half of the chain that 
had an even number of sites. 
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FIG. 11. (Color online.) The wave number g ost against J2. 
Data is shown for four chains of different length (a) and with 
the result of a fit q CEt = 1.2062 9(J 2 - 0.528)(J 2 - 0.528) ' 4806 
for a chain of 601 sites (b). 



To estimate the wave-length of the incommensurate 
modulation we make use of the fact that, if our system had 
translational invariance, the distance between maxima 
in the on-site magnetization would be equal to half of 
the wave-length, as indicated in Fig. [To[e). Thus, by 
calculating the mean distance of the central maxima, we 
are able to determine an estimate for the wave-length 
of the incommensurate modulation. The inverse of this 
quantity can then be used to calculate the wave number, 
2n/X cst . In Fig. 11 a) we show how q ost varies with 



sites. 

Since the incommensurate behavior can only be seen 
if the wave-length is shorter than the system, it starts 
later in smaller chains. Aside from small deviations, that 
can be attributed to finite size effects, the wave-length 
does only depend on J 2 and not the length of the chain 
(see Fig. [TT|a)). In the limit of infinite J2, the next- 
neighbor interaction J\ can be neglected and the chain 
be partitioned into two sub-chains that do not interact. 
As mentioned, the Jy J2-model in this limit approaches 
two uncoupled chains with intra chain coupling J2. The 
wave-length of the incommensurate behavior in this limit 
reaches its minimum with A = 4 lattice spacings. 

For J 2 > J 2 L , one expects the wave number q to be- 
have as q oc (J2 — J^) : j where < a < 1.P In a study 
of correlations functions around the disorder point in 
the S — 1/2 Ji-J^model with an even number of sites 
and modified interactions on the edge of the chain, the 
exponent was reported to have been calculated to be 
a = 1/2 JS2l This i s con sistent with calculations on classi- 
cal Lifshitz points^J^SI which at the mean field level find 
a = 1/2. In the present case, where the ground state 
momentum is changing, one might also expect corrections 
to the mean- field value of a = 1/2 as described in Ref. HT1 
and SSI 

Our calculations indeed confirm that q es t(J2) follows 
a power-law with exponent smaller than 1 (see Fig. [TT| 
b)). The line in Fig. 11 'b) is a fit of the three parameter 
function f(x) = ci<3(x — c%)\x — c 2 |", where Q{x) is 
the Heaviside step function, to the blue data points also 
shown in the plot. Using this form we find a value for the 
exponent a — 0.4806. The data in Fig. [XT] show steplike 
features. The cause of the steps is the introduction of new 
maxima: every time a new maximum appears, g e st jumps 
abruptly in order to accommodate the new maximum and 
there is a step. Between the appearance of new maxima, 
the maxima that are present move closer together and g cst 
increases smoothly. As one increases the system size, this 
effect affects the mean distance between maxima less and 
thus leads to less pronounced steps. Due to the different 
range of ^-values, the steplike features explained above 
are more pronounced in Fig. [TT|b) than in Fig. 11 a). 
Because of the inaccuracies the steplike features introduce 
to the fitting procedure, we cannot comment on whether 
or not the corrections mentioned above are necessary. 

While the on-site magnetization could relatively easily 
be understood from the results obtained with PBC, this 
is not the case for the energy spectrum. To the left of 
the transition (J 2 < J^), the spectrum for OBC (shown 



lest 

J2 for four chains whose length ranges from 301 to 1001 



Fig. 12 1 looks exactly like the spectrum for PBC (shown 
Fig.jTb - yet there is an important difference: the spectrum 
for OBC is not degenerate. Introducing the boundary 
splits the degenerate states. On the other side of the 
transition (J2 > J2), the behavior of the energy of the 
first excited state also looks familiar: it hits the ground 
state energy, grows, approaches it again and another 
level-crossing occurs. Repeated level-crossings of just the 
two states follow. Higher excitation- levels, however, do 
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defined point in the spectrum. In order to answer this 
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FIG. 12. (Color online.) The excitation spectrum with open 
boundary conditions in the variational subspace with Sj ot = 
1/2. The energies of the first few excited states are shown. 
The ground state energy was set to zero. The values at which 
the scaling with N is studied in Fig.[l4]are indicated by dashed 
vertical lines. The data was taken for a chain with 601 sites. 



not cross many other levels as they do for PBC. They 
approach the ground state, then turn around and form a 
pair with the state they would have been degenerate with 
under PBC. The two states exhibit a repeated pattern 
of intertwining level-crossings while their mean energy- 
difference to the ground state grows. We do not know of 
an intuitive way of understanding the spectrum for OBC 
from the spectrum with PBC. Modifying the couplings 
at the boundary of the chain by a multiplicative factor 
of A and varying A between and 1 we have studied the 
cross-over from PBC to OBC. A low-energy spectrum 
similar to the one for OBC is observed until A w 0.9. 



question we show results in Fig. 13 for chains of length 
301 and 451 sites for a range of J 2 close to . As can be 
clearly seen, the minima of the higher energy levels occur 
much closer to the first level crossing of the ground state 
for N — 451 than for N = 301. In the thermodynamic 
limit we expect the minima for all higher lying levels to 
occur at 

We next focus on the scaling of the energy levels with N. 
In Fig. 14 'a) we show data for N 2 E n taken at J 2 = 0.526, 
to the left of the transition as indicated in Fig. [12] and [13] 
As can be seen, it converges to a constant value indicating 
that for this value of J 2 E n cx A -2 . For the first excited 
state this behavior is apparent for quite short chains al- 
ready and it seems plausible that for higher excited states 
longer chains would lead to the same decay proportional 
to N~ 2 . This scaling is not surprising since the soliton 
behaves like a massive particle in a box. We therefore 
expect the low energy spectrum to be approximated by 
h 2 k 2 /(2A sol ) with k = nn/2N, n = 1,2.., yielding the 
expected scaling of the energies as A~ 2 . 
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FIG. 13. (Color online.) Spectrum with open boundary condi- 
tions for chains of 301 and 451 sites. The values at which the 
scaling with N is studied in Fig. [14] are indicated by dashed 
vertical lines. 



It is reasonable to ask if the Lifshitz point is a well 



FIG. 14. (Color online.) The energy scales proportionally to 
N~ 2 . Plot of N 2 AE n to the left of the transition (a) and to 
the right of the transition (b) for the first five excited states. 

In Fig. [l4^b) we show data taken on the other side of 



J ^ at J; 



12 



and 



13 I. For 



2 = 0.531 (again indicated in Fig. 
the smallest N shown, the higher excited states still show 
signs of the transition at this value of J 2 . For short chains 
the second, third, forth and fifth excited state thus have 



minimum in Fig. 14 b). For chains with more than roughly 
160 sites we see of intertwining pairs of states familiar 
from Figs. |12| and |13| The average energy of the pair at 
big N also scales proportionally to N~ 2 . We therefore 
conclude that sufficiently far away from the transition 
point for the first few energy-states AE n cx J 2 /A 2 . 

In order to study the scaling of the spectrum at the 
Lifshitz point we focus on the minimum in the sec- 
ond excited state. Although this minimum occurs at 
slightly different J 2 as N is varied it serves as the best 
possible definition of an excited energy scale at the Lif- 
shitz point. Specifically, we define the minimal energy- 
difference of the ground state and the second excited 
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FIG. 15. (Color online.) The energy at which the 2 nd excited 
state turns around goes to zero faster than N~ 2 . The value 
of J 2 of E? in was found up to AJ 2 = 10~ 5 . The resulting 
uncertainty of the value of the minimum energy is smaller 
than the size of the symbols in the plot. 



state as Ef m 



Our results for Ef m 



min, / 2 [E 2 - E 
are shown in Fig |15| As can be clearly seen in this figure 
£™ ln goes to zero faster than N~ 2 violating the simple 
scaling found elsewhere. 
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FIG. 16. (Color online.) Main panel: J% as a function of the 
length of the chain. Inset: The difference between J 2 and the 
C2 at which the minimum of the second excited state occurs. 
All values were determined up to AJ2 = 1CP 5 . The resulting 
uncertainty of the value of the minimum energy is smaller 
than the size of the symbols in the plot and causes deviations 
from very smooth behavior. 

We now turn to an estimate of the location of J^ within 
the variational approach. The level-crossing of the first 
excited and the ground state allows for an easy way to 
define the value of Jk for a given JV. As one could already 

J2 varies slightly with the 



see in the Figures 12 and 
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length of the chain. Our results are shown in Fig. [Id] for 
chains out to JV = 701. The main panel in Fig. [T6| shows 
that J2 converges to approximately = 0.528 as one 



increases the length of the chain. 

The value of J2 at which the second excited state has 
its minimum also approaches J% . To show this we use 
the value of J2 at which the n-th state reaches its first 
minimum fo r a given JV. We call this quantity C n . The 
shows AC21 = C2 



inset in Fig. 



16 



J 2 ■ 



As can be seen, 
this quantity approaches and the minimum for big JV 
thus lies at the Lifshitz point. 



IV. 



INCOMMENSURATE BEHAVIOR IN THE 
ENTANGLEMENT ENTROPY 
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FIG. 17. (Color online.) The entanglement entropy of a 
bipartition of an odd chain. The data was taken for chain of 
301 sites and the variational subspace Rq (see Sec. [TTJ> was 
used. 

The scaling of the entanglement entropy at a (quan- 
tum) Lifs hitz t ransition has recently been the subject 
of interest!- 14 * 15 ! In free fermion models, analogous to the 
spin chain model discussed here, the Lifshitz transition 
is associated with a change in the topology of the Fermi 
surface. In one dimension new Fermi points appear at the 
Lifshitz transition and, analogously, new patches appear 
in higher dimensional models. If one associates a chiral 
conformal field theory with each patch, it can be arguecP^ 
that, when the number of points (patches) increases by 
a factor K, the entanglement should be multiplied with 
the same factor K . For a free fermion model with next 
nearest neighbor hopping, £2, one expects the number of 
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FIG. 18. (Color online.) The entanglement entropy of a 
bipartition of an even chain (200 sites) (a) and an odd chain 
(201 sites) (b). The data was obtained using DMRG with 
m = 256 states kept. 



Fermi points to double at the Lifshitz transition t% = 1/2 
at half-filling with a corresponding doubling in the en- 
tanglement entropy. This behavior is well confirmed in 
numerical calculations.^ 

In this section we discuss our results for the entangle- 
ment entropy across the Lifshitz point in the odd length 
J\-Ji quantum spin chain which is the quantum spin 
analogue of the model considered in Ref. [T3J We study 
the entanglement in terms of the von Neumann entan- 
glement entropy of a sub-sys tem A of size I and reduced 
density-matrix p A defined bjEMl, 



S(l,N) = -Ti[p A lnp A ] 



(12) 



where N again stands for the total system size. We 
consider exclusively open boundary conditions. 

If one uses the restricted space Rg, which was intro- 
duced in Sec. [TTJ, as the variational subspace, one can 
also calculate the entanglement entropy using the method 
employed in this paper psl Away from MG and Lifshitz 
points, where the variational method is not reliable, we 
complement the variational results with DMRG calcula- 



tions. 

We first discuss the variational results for odd length 
chains close to the Lifshitz point shown in Fig. [17] for 
N = 301. The entanglement entropy at the MG point for 
the odd length chain, shown in Fig. 17 a), has previously 
been discussed in detail. 28 Since the entanglement en- 
tropy is very directly connected to the wave function of 
the state, drastic changes of the wave function should also 
be present in the entanglement entropy when the Lifshitz 
poin t is reached. This is clearly the case as can be seen in 
Fig. 17 As the Lifshitz point, J^, is reached, the entan- 
glement entropy develops plateaus (Fig|17[c)). As J2 is 
increased more plateaus appear (Fig. 17 d)j. For the free 



fermion model studied in Ref. [To] analogous oscillations in 
the entanglement entropy are observed beyond > 1/2. 
Because a different subspace was used in the previous 
parts of this paper, the transition begins at J 2 ~ 0.529 
which is slightly higher than J2 ~ 0.528 which could be 
inferred from Fig. [13] 

For an even length system no such plateaus are visible 
(see Fig. |18[ a)). As J2 — > 00 the entanglement increases 
towards that of two independent gapless Heisenberg chains 
as it must. A similar increase is seen for an odd number 
of sites but with pronounced signatures of the incommen- 
surability (see Fig. [l8^b)). 



V. THE TRANSITION POINT 

The numerical value of J 2 for the Lifshitz point depends 
not only on the length of the chain but also on the basis- 
set that one uses in the variational calculation. While this 
is a small concern when one looks at qualitative features, 
it is of course detrimental if one is interested in a precise 
estimate of the Lifshitz point. Just using the different 
basis sets introduced in Sec. [TT] this is evident. Using 
the smallest basis, Rq, for a chain with N = 301 sites, 
we obtained 



0.5295 (see Fig. [17| for the onset of 
oscillations in the entanglement. This is a slightly bigger 
value than what was found in Fig. O J 2 L w 0.528, based 



on the calculations with the larger basis Z s . 

DMRG can give us a more reliable estimate for at 
least for small chains. For a chain with 201 sites we found 
the first indications of incommensurate behavior in the 
local magnetization at: 



0.538(1) 



(13) 



We expect this estimate to depend on N in roughly the 
same way as the variational estimate does in Fig. [16] If 
this is the case, an eventual extrapolation to the N — > 00 
limit might change this estimate by 0.0005 which is smaller 
than the uncertainty to which we have determined the 
point. 
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VI. CONCLUSION 

We have studied incommensurability effects as they 
occur in the odd length antiferromagnetic J1-J2 chain. 
Even though no exact ground state wave function is known 
at the MG point, J2 = Ji/2, this point is the disorder 
point with minimal correlation length. The Lifshitz point 
J% = 0.538Ji marks the onset of significant modulations 
directly in the ground state (Sf) as well as a shift in the 
ground state momentum. A series of inter-twining level 
crossings causing the shift in the ground state momentum 
starts at the Lifshitz point. The shift in the ground state 
momentum and the associated modulations directly affect 
the entanglement entropy which shows distinct plateaus 
developing for J 2 > J%. 

In realistic compounds with chain breaking impuri- 
ties one would expect half the chain segments to be of 
odd length. The experimentally well studied compound 



CuGe0 3 has aJ 2 ~ 0.36 Ji < J1/2P If compounds with 
a J2 in excess of Ji/2 can be identified, it would be very 
interesting to experimentally look for the odd length ef- 
fects that we have detailed here. In particular, the effects 
on the on-site magnetization shown in Fig. |10| might be 
observable using NMR techniques or other local probes. 
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